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Here we apply our SU(N) and U (N) parameterizations 0, Q to the question of entanglement 
in the two qubit and qubit/qutrit system. In particular, the group operations which entangle 
a two qubit pure state will be given, as well as the corresponding manifold that the operations 
parameterize. We also give the volume of this manifold, as well as the hypothesized volume for the 
set of all entangled two qubit pure and mixed states. Extension of this work to the qubit/qutrit 
system will also be given. 
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I. INTRODUCTION 

Wc know that, in general, TV-dimensional density matrices can be written in the form 



P=4 W^V*)]-' (D 



N\ V 2 
This representation is a convenient one since for pure states it yields, 

n ■ n = 1, and n * n = n, (2) 
where the "star" product is defined according to 



(a*b) fe = \j N N ^ ( jdijkdibj. (3) 



2 \N -2 

The way to see this is to just use the following relationship from 0, [j| 

2 

XiXj = j^S^In + (dijk + icijk)\k, (4) 

where the Cijk and the dijk are the structure coefficients for the SU(N) Lie algebra in question. For TV-dimensional 
density matrices there exists a simple procedure to calculate the components of n: beginning with equation we 
use the trace condition on the elements of the algebra 

Tr[\i • A,] = 2% (5) 

to generate 



2N{N- 1) 



* Tr[(Np — t N ) ■ Xi 



A *Tr[p- Xi] (6) 



2(N- 1) 



We can use this to represent n in terms of other parameterizations such as when p = U pdU^ 2 or to evaluate Q when 
p is explicitly given as in the following example. 

A. Example Calculation: Bell's States 

We begin by defining the four Bell states @ : 






= -?= n h ^2 = — = n , ^3 = -7= 1 , V^4 = i I 3 (7) 




1 Note that orthogonal states in this representation are those that have an angle 9 = cos 1 ( ) between them in the corresponding 
Hilbert space. 

2 Substitution of p = Up^U' 1 into equation JfjJ yields 



TV * 

Tr[(t/p d i7t) • Xi] 



2(N - 1) 



which upon calculation gives n as a function of the TV (TV — 1) group parameters (denoted by a,;) and the TV — 1 "rotations" (denoted by 
8i) parameterizing p^. 
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which can also be represented as the following density matrices 



4>ii>\ = pbsi 



■03 ^3 



PBS3 



( 1 1 





V i o o i 

/ 

110 

110 

V o o o o 



1p 2 ^l = PBS2 



PBSi 



/ 

I 


1 
1 


n 
u 


n 

u — 


U 


U 


u 













\ 


-1 




















1 




-1 


( 


- 


1 


1 













(8) 



These four density matrices represent the four possible EPR pairs for two qubit systems in SU (4) and represent an 
orthonormal basis for the entire two qubit, pure state, state space (i. e. for the vector space but not for p). By their 
definition, they are maximally entangled states, i.e. nonfactorizable superpositions of product states, and thus impart 
non-local correlations between the behavior of the two qubits that make up these four states. For, if these states were 
factorizable 



P = PA® PB, 

then the probability of joint detection would also factorize 

\p\ 2 = \pa\ 2 \pb\ 2 , 



(9) 



(10) 



and thus the measurements would be independent of each other. 

Now although the Bell states cannot be decomposed into a set of product states, we can decompose their density 
matrix representation given in equation (JSJ into the one given by equation (JTJ by using equation JBJ) with N = 4, 



Vein - a), 



(11) 



thus yielding 



Pbsi 



PBS2 



PBS3 = 



PBSi 







3V2 A8 + 




^A 15 )), 




,V6(-La 3 + 


3V2 XS 


^- 


^A 15 )), 






+ y|Aa- 


3V2 A8 


+ ^A 15 )), 






-y|A 6 - 


3V2 A8 


+ ^A 15 )). 



(12) 



We can see that for the first two Bell states the only non-zero components of n are n^, ns, ng and nis whereas 
for the last two Bell states one switches ng for ne (with its corresponding Lie algebra component). The rest of the 
components of n are zero. 

II. GENERAL ENTANGLING OPERATIONS ON TWO QUBITS 

From (OiBEJ and references within) we know that entangled (and entangleable) two qubit density matrices must 
satisfy 3 



Tr[p 2 } > o- 



(13) 



Note that this is just a necessary condition for a state to either be already entangled, or to possibly be entangled under some group 
operation. Obviously, pure states satisfy this criterion automatically, but it is the mixed state situation to which this criterion is more 
often applied. Separable two qubit density matrices have Tr[p 2 ] < ^. 
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If we start with equation , assigning N = 4 



\ (li + V6(n • A) 



then equation (|13|1 yields 



Tr[p 2 } = Tr 



^ (i 4 + V6(n ■ a; 



2-i 



^(l + 3||n|| 2 )>i 



(14) 



(15) 



implying that ||n|| > 1/3 for entangled states. On the other hand if we begin with p = UpdU^ where U £ 5*J7(4), we 
can see that 



Tr[p 2 ] = Tr[(Up d tf)(Up d tf)] = Tr[Up 2 U^} = Tr[p 2 d ] > ~ 



where pd is 



Pd 



/sin 2 (6>i)sin 2 (6l 2 )sin 2 ((9 3 ) 

cos 2 (6>i)sin 2 (6>2)sin 2 (6> 3 ) 



V 






cos 2 (0 2 ) sin 2 




\ 




cos 2 (0 3 )/ 



(16) 



(17) 



from [lj and we have exploited the knowledge that UU* — I4. Evaluation of this trace yields the following demands 
on the ranges of 9\, 2 , and 63 (here expressed as functions of sin(#i) and cos(6>i)): 



Tr[p 2 d ] = cos(0 3 ) 4 + (cos(0 2 ) 4 + 1 (3 + cos(4 0{)) sin(# 2 ) 4 ^ sin(0 3 ) 4 > I. 
Depending on the ranges of 9i one could either have a pure or mixed state that would satisfy equation l|13fl • 



(18) 



A. Pure State Entanglement 



One can see that for Oi = tt/2, (i = 1,2,3), equation (|18|l would be satisfied and thus we would have a pd which 
could be entangled. Therefore using equation H17|l with 9{ = 7r/2, (i = 1, 2, 3), we can define a pure state as 4 



Pd 



(\ 0\ 





Vo 0/ 



(19) 



Using equation © we can then calculate the components of n in terms of the twelve a and three 9 parameters when 
p is given by p = UpdU^ and U € S'C/(4). From this it would be possible to determine the actual parameter values 
that would generate the Bell states given in equation © by solving the fifteen simultaneous equations implied by the 
representation of n in terms of the Euler parameters. Another, more instructive way is to apply successive unitary 
operations, U € 5J7(4), to the two qubit pure state given in equation i|T§|) until one achieves the requisite Bell state. 5 



4 Remember that any 4 by 4 matrix with one element along the diagonal equal to unity and the rest equal to zero will be invariant under 
a U(3) subgroup of SU(A) and would therefore represent a pure state. The representation given in equation 1191 however is ideally 
suited for the following calculations based on the group parameters remaining after evaluating the coset SU(A) /U (3) (see Q for more 
details) since it is those remaining 6 parameters (equally split between three A3 phases and three Aj i _ 1 j2 +1 (i = 2,3,4) rotations) 
that are explicitly contai ned in the pure state volume measure given in which parameterizes the space CP 3 via the corresponding 
Fubini-Study metric (see lilt ll ll for the SC/(3)case) which directly acts upon the P-f^i} element. This pure state representation is also 
consistent with the generalized p^i used in [UlJ given previously. 

5 Because of the U(3) invariance of the two qubit pure state given in equation 1191 . we only have to look at those operations in the coset 
SU(A)/U(3) which we know from to be represented by 

gS^3 Q l ^^2^2 e i^3 a 3 e i ^5 a i e i ^3 a 5 gi^W^B 

but it is instructive to see, through the use of the full group SU (4), that only those operations generated from the set {A2, A5, A10} yield 
pure state entanglement. 
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Bell States One and Two 



To begin, we first act upon our pure state with the group operations, 



U 



oa and £7 f 



-i\ 10 a 



yielding 



which in matrix notation is 



p = U Pd U^ = e iX ™ a p d e 



cos(a) sin(a) N 

10 

10 

- sin(a) cos(a) 

/ cos(a) 2 



/l ON 





\0 0/ 



iAioc* 



/cos(a) 
1 
1 



\sin(a) cos(a) 




cos(a) sin(a)\ 




\— cos(a) sin(a) sin(a) 2 / 



Taking the partial transpose of the above density matrix yields 



/cos(a) 2 



pOt = 



V 





— cos(a) sin(a) 
cos(a) sin(a) 

sin(a) 2 ■ 



which has an eigenvalue decomposition equal to 

{Xi,X2,X3,Xi} = {cos(a) 2 , -cos(a)sin(a), cos(a) sin(a), sin(a) 2 }, 
and where the constant term in the characteristic polynomial is 6 

Bct(p pt - I4 * x) ->■ - cos(a) 4 sin(a) 4 . 



(20) 
(21) 



(22) 



(23) 



(24) 



(25) 



Recalling that a varies from to tt/2 we can see that for < a < % we have an entangled density matrix p. In 
particular we can see that if a = ir/4 we generate the second Bell state 



cos(f)sin(f)N 



sin(f) 2 



/ cos(f f 

tt\ I 

AJ~ 

V-cos(f)sin(|) 

(\ -ft 



V-3 3 J 

Now, in general, if we use the Euler angle parameterization of SU(4) given in Q and take the 
U £ SU(4) to be given as 

U = U[cc\, 0, «3, 0, a 5 , a 6 , ar, 0, ctg, 0, an, 0, 013, au, a 15 ) 

_ e iA 3 ai e iA3Q3 e iA3Q5 e iAioQ;6 e iA3Q7 e iA3ag e iA3aii e iA3ai3 e iA8Q!i4giAi5ai5 



(26) 

most general 
(27) 



The constant term is just the zeroth order coefficient of \ in the characteristic equation. 
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we would then have a density matrix equal to 
P = Up d U^ 



t 



cos(a 6 ) 2 -e l ( Ql + Q3+Q5 )cos(a 6 )sin(a 6 )\ 





< ai+a3+a ^ cos(a 6 )sm(a 6 ) sin(a 6 ) 2 



(28) 



whose partial transpose is 



/cos(a 6 ) 2 

-e 

-e i ( ai + Q3+Q5 )cos(a 6 )sin(a 6 ) 
\ 



\ 

i(ai+Q3+a5) cos(a6)sin(a 6 ) 



sin(a 6 ) 2 / 



(29) 



which yields an eigenvalue decomposition and a constant term in the characteristic polynomial equivalent to equations 
(|24|l and 1)25(1 but that does not generate the second Bell state when a§ = it/ 4 



cos(f) 2 





-e^ Ql+a3+Q5 )cos(f)sin(f)\ 








. e -i(«l+O3+O6) c0s (|) 8 i n (|) 






sin(f) 2 



-l e ^ ai+a3+a ^\ 





(30) 



/ 



What we would have, on the other hand, is the following. Define ct\ + + a§ = j3, where < (3 < 57r since ct\ runs 
from to 7r and «3 and «5 each run from to 2tt7 Then 

e ±i(«l+Q3 + Q5) _ g±«/3 

= cos(/3) ± i sin(/3) 

= — 1 when P = 7r, 37r, 5n 

=1 when (3 = 0,2vr,47r. (31) 

Thus when ai + a 3 + a 5 = n, 37r, or 5ir we get the first Bell state and when u\ + + a$ — 0, 27r, or Air we get 
the second Bell state. Intermediate values of u\ + + as can be equated to "intermediate" Bell states; states which 
have an equivalent density matrix representation as the first and second Bell state, but which are not equal to any 
type of convex sum of said Bell states. 

Now if, instead of equation i(2U|l we were to choose 

U = e zX2a orU= e lXsa (32) 

and apply it to pa (as in equation (|21|) ') we would not generate the other two Bell states, or for that matter any entangled 
density matrix p, even with the most general U (as in equation (|27(l 'l. The question now is, what combination of the 
exponentiation of A2 with A5, and/or A10 will entangle the pure state density matrix p^ yielding the other two Bell 
states. It is to this question we now proceed. 



2. Bell States Three and Four 



To begin, we first act upon our pure state with the group operations, 

U = e lX ^e lX2U and = e^^e"^, (33) 



' We are now using the covering ranges for 5(7(4), denned in Q. 
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yielding 



U Pd U^ 



(34) 



which in matrix notation is 8 




























Hi 








°, 


( 










0/ 





cos(/z) cos(^) cos(/x) sin(i/) sin(/i) 

— sin(z/) cos(^) 

— cos(^) sin(/i) — sin(/x) sm(u) cos(p) 
1 

/ cos(p,) 2 cos(f ) 2 — cos(/x) cos(i/) sin(i/) — cos(/x) cos(;/) 2 sin(/x) 

— cos(/i) cos(^) sin(z/) sin(y) 2 cos(v) sin(/z) sin(z/) 

— cos(p) cos(i^) 2 sin(/x) cos(i^) sin(//) sin(^) cos(z^) 2 sin(^) 2 
V 



cos(^i) cos(^) sm(v) cos(v) sin(/z) 0\ 

- cos(n) sin(i/) cos(^) — sin(/i) sin(z^) 

— sin(/x) cos(p) 





(35) 



Now the above density matrix does not look like either of the remaining two Bell states; unless we demand that 
fi = n/2. Then we get 



'0 0\ 

sin(^) 2 cos(^) sin(^) 

cos(^) sm(i>) cos(i^) 2 

v0 0/ 



(36) 



which has the same form as the remaining two Bell states. As before, taking the partial transpose of the above density 
matrix yields 



/ cos(^) sin(i/)N 

sin(^) 2 ' 

cos(^) 2 

\cos(z/) sin(z/) 



which has an eigenvalue decomposition equal to 

{xi,X2,X3,X4> = {cos(i/) 2 , - cos(v) sin(i/), cos(z/) sm(v), sin(j/) 2 }, 
and where the constant term in the characteristic polynomial is 

Dct(p pt - 1 4 * x) -> - cos(^) 4 sin(^) 4 . 



(37) 

(38) 
(39) 



Recalling that v varies from to 7r/2 we can see that for < v < ^ we have an entangled density matrix p. In 
particular we can see that if v = n/4 we generate the third Bell state 



P \ v = T 



/0 0\ 

sin(f) 2 cos(f)sin(f) 
cos(f)sin(f) cos(f) 2 



V0 











0/ 



Notice, if we had instead used 



/0 0\ 

= o \ \ o 

I I ' 

\0 oj 



(40) 



(41) 



Recall that since A2 and A5 do not commute, the other possible group operation e lA2 "e lAs,l p £ je ,A 5M e »A 2 f w ;n no t generate this 
matrix. But, because [A2,As] = — [Ab,A2], the subsequent work after this step will be similar for either e lA2I/ e lABM pde~ lA5M e~ lA2 ' y or 
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as our initial starting point, we would have instead produced the following density matrix 

cos(cj)) 2 cos(V') 2 — cos(0) 2 cos(V') sin(^) — cos(0) cos(-i/>) sin(^)) 0\ 

— cos(</f>) 2 cos(i/0 sin(-0) cos(0) 2 sin(-0) 2 cos(0) sin(0) sm(ip) 

— cos((/>) cos(V') sin(0) cos(0) sin(0) sin(V>) sin((/>) 2 
0/ 

One can see that in order for this matrix to be equivalent to l|36|) . one must demand that ip — 7r/2, thus yielding 



(42) 



/0 0\ 

cos(cf)) 2 cos(^) sin(0) 

cos(^) sin((/)) sin(0) 2 

Vo 1 o o oy 



(43) 



For simplicity, we shall generalize this group operation and not the original e^^e 1, 2V calculation. 

Therefore, in general, if we use the Euler angle parameterization of <SZ7(4) given in p| and take U £ SU(4) to be 
given as 



U = U(ai, -,a 3 ,Q!4,Q;5,0, a 7 ,0, a 9 ,0, an,0, ai 3 ,ai4,ai 5 ) 

__ e i^3ai e i\2^ e i^3^3 ^i^^ai^iXaa^ ^i\ 3 a 7 ^i\ 3 a a ^i\ 3 aii ^1X30113 ^iXgati^iXt^a! 



(44) 



we would generate the following density matrix 

P = Up d rf 





cos(a4) 2 e" 

e l ( Ql - Q3 )cos(a 4 )sin(a 4 ) 

\0 



0\ 

-»(ai-a 3 ) cos (q; 4 ) sin(o;4) 

sin(a 4 ) 2 

0/ 



(45) 



whose partial transpose is 
/ 



p 



e l ( ai - Q3 )cos(a 4 )sin(a 4 )^ 

cos(a 4 ) 2 

sin(a 4 ) 2 

e-^ ai - a3 )cos(a 4 )sin(o!4) / 



(46) 



which yields an eigenvalue decomposition and a constant term in the characteristic polynomial equivalent to equations 
(|38ll and l|39|) but that does not generate the third Bell state when a 4 = 7r/4 



P\ «4 



A) 












o\ 





cos(f) 2 






-<*d cos(f 


)sin(f) 





e< ai - a3 )cos(f)sin(f) 






sin(f) 2 





\o 














fo 





°^ 











1 Ip— »(ai— as) 

2 2 














l p i(ai-a 3 ) 1 
2 C 2 











v° 














What we would have, as before with the second Bell state, is the following. Define a% 
since ol\ and 013 again run from to 7r and to 2-7T respectively. Then 



(47) 

a 3 = 7j where —2tt < 7 < 



e ±i(a 1 ~a 3 ) _ e ±i7 



= cos(— I7I) ± i sin(— I7I) 
= cos(|7|) =F isin(|-y|) 
= — 1 when I7I = 7r 
=1 when I7I = 0,27r. 



(48) 



Thus when \a± — = tt we get the fourth Bell state and when |ai — as| = or 2ir we get the third Bell state. 
Intermediate values of \a.\ — can be thought of as "intermediate" Bell states; states which have an equivalent 
density matrix representation as the third and fourth Bell state, but which are not equal to any type of convex sum 
of said Bell states. 
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3. General Two Qubit Pure State Entanglement 

The natural extension of the previous work is to look at the case when we use the Euler angle parameterization of 
SU(4) given in [IJ and take U £ SU(4) to be given as 9 



U = U (ot x , a 2 , a 3 , a 4 , a 5 , a 6 , a 7 , 0, a 9 , 0, an, 0, a 13) a 14 , a 15 ) 

_ giA3Qi giA2«2 giA3Q3 ^1X5014 giA3Q5giAioa6giA3a7giA3a9giA3Q;ii giA3ai3 ^lAga^ giAisais 



(49) 



We would then generate the following density matrix (see p on next page). One can see immediately that in order to 
obtain the general form for Bell states three and four, «2 and must be set to 7r/2 and zero respectively. Similarly, 
in order to obtain the general form for Bell states one and two, ui and ay must be set to zero. 
Generally, though, the two sets of eigenvalues of the partial transpose of equation (|50jl 

*± = ± e ^( 2Q i+^) cos ( a4 ) CO s(a 6 ) x A, 



1 ± i e T i (2ai+a 5 ) / ^ e i(2Q 1+ a 5 ) _ 4cos(a 4 ) 2 cos(a 6 )) X A, (51) 
where A is 

i 

A = (e M5 cos(ae) sin(a2) sin(a4) + e 2lai cos(«2) sin(ae)) 2 

i 

x (e 2lQ1 cos(a 6 ) sin(a 2 ) sin(a 4 ) + e la5 cos(a 2 ) sin(a 6 )) 2 , (52) 

indicate that the phase parameter 03 does not contribute. We claim that one only needs three rotations, and one 
overall phase in order to carry out a general entangling operation on pd- This can be seen if one expands the constant 
term from the characteristic polynomial for this situation 



Det(p J,t - I4 * e) -»• - e - 2 *( 2 «i+«5) cos ( a4 )4 cos ( af .) 

x (e IQs cos(a 6 ) sin(a 2 ) sin(a 4 ) + e 2lQl cos(a 2 ) sin(a 6 )) 



4 

2ia, / s . / n\2 



x (e 2M1 cos(ae) sin(a2) sin(a 4 ) + e las cos(a2) sin(ae)) 2 (53) 



and defines r\ = 2ai — as as the cumulative phase 



e° = — cos(a 4 ) 4 cos(a6) 8 sin(a2) 4 sin(a 4 ) 4 

— 2e~ l?, cos(a2)cos(a 4 ) 4 cos(a6) sin(a2) 3 sin(a 4 ) 3 sin(ae) 

— 2e" 7 cos(a2)cos(a 4 ) cos(ag) sin(a2) 3 sin(a 4 ) sin(ag) 

— 4cos(a2) 2 cos(a 4 ) 4 cos(a6) 6 sin(a2) 2 sin(a 4 ) 2 sin(a6) 2 



2 2lr? cos(a2) 2 cos(a 4 ) cos(a6) 6 sin(a2) 2 sin(a 4 ) 2 sin(a6) 2 



e 2l?J cos(a2) cos(a 4 ) cos(ae) sin(a2) sin(a 4 ) sin(ae) 
e 2 

— 2e _ "'cos(a2) 3 cos(a 4 ) cos(ae) 5 sin(a2) sin(a 4 )sin(ag) 3 

— 2e lr? cos(a2) 3 cos(a 4 ) 4 cos(a6) 5 sin(a2) sin(a 4 )sin(ae) 3 

— cos(a 2 ) 4 cos(a 4 ) 4 cos(a 6 ) 4 sm(a 6 ) 4 . (54) 

From the previous work it is obvious that only the parameters a2i (i = 1,2,3) and 77, the overall phase running 
from to 2ir, are needed to parameterize all entangling operations that can be done on an initial pure state (given in 
equation 119(1 ). These three rotations and one phase thus can be used to parameterize the following manifold 



51/(4) _ CP 3 

CT(3) x U(l) SU (2 )n ~ U(l) su(2)2ai _ a5 



M e S[/(4) = TT(0 ^ z TT7T\ = Tr7r\ > ( 55 ) 



9 We only need to look at the first 6 group operations of U because they are the ones that "parameterize" the coset SU(4)/U(3) = CP 3 



p = Up d U x = 

cos(ci2 J 2 cos ((14) 2 cos(ojq ) 2 — e 2 ^ a l cos(a2 J cos(ct4) 2 cos(ojq ) 2 sin(ct2 J e**- a l~^ a 3) cos (a 2 ) cos(cKg ) 2 sin ( 2a 4 ) — —J- e i *- cv l~'~ a 3~'~ a 5) cos(a2 ) cos(ci4 J sin(2ag ) 

— 2ia 1 cos ^q, 2 ) cos(ck4 ) 2 cos(ckq) 2 sin(ct2 ) cosfa^) 2 cos(ag) 2 sin (a 2 ) 2 e — ~~ a 3 ) cos(c*4) cos(ckq) 2 sin(a2 ) sin (0^4) e~ *( a l ~~ a 3 — a 5 ) cos (a 4 ) cos(ag ) sin (a 2) sin(aQ ) 

H Q 1+ Q 3) cos(q2) cos(ct4) cos(ojq) 2 sin(ci4) e^^ a ^ cos (£14) cos(ckq) 2 sin (0:2) sin (a 4 J cos(cKg) 2 sin (0C4) 2 e ia 5 cos(ag) sin (04) sinfojg) 

o 1 +a 3 +a5) cos(a2) cos(a4) cosfcKg) sin(ag) e*' * 1 — a 3 — a 5 ) cos (04) cos(ckq) sinfa^) sin(ae J e~~ ia 5 cos(qq) sin(a4) sin(ag) sin(cte J 2 



(50) 



11 



the volume of which can be found using the material from 



wt4 » v^(3) xc/(i) S£/ (2),y v^(i)2*a 3 ; 2 * v^ cl)X3 

_ 7T 3 /6 

~~ 2?r x 2tt 

= (56) 

What we have been able to derive in equation (|56(l is the volume of the manifold of all operations on two qubit pure 
states which produce entanglement. Since these operations act upon a pure state (again, see equation lf[T?|> ). which is 
just a point in CP 3 and thus of measure zero, one may conclude that, up to the volume of a measure zero set, the 
volume of the manifold of all operations on two qubit pure states which produce entanglement is equivalent to the 
volume of the set of all entangled two qubit pure states 



{Set of All Entangled Two Qubit Pure States} = Me su 



(4) 



7T 

Vket of All Entangled Two Qubit Pure States = ^Me s „ (4) = ^J- (57) 

This volume is less than the numerically estimated value calculated by Zyczkowski et. al. in Q(eq. 30) and referenced 
in for the lower bound of the volume of entangled two qubit states (1-0.863) by approximately 6 one-thousandths. 

Unfortunately, whereas our calculation was done with only pure states in mind, the Zyczkowski et. al. calculation was 
numerically done using a mixed state product measure defined on U (AT) and randomly chosen density matrices H, 
which satisfied the Peres-Horodecki criterion for separability 0, 0] . Therefore although a tantalizing conclusion, we 
must concede that until a more general calculation is performed using our mixed state product measure (defined in 
0) we cannot definitively state that we have calculated the exact volume of the set of all entangled two qubit pure 
states. 



B. Mixed State Entanglement 



Using common 
satisfied if 



or 



or 



mathematical software, one can see that for di ^ 7r/2, (i — 1,2,3) equation (|18|l could still be 



. 1 1 V3278975 = 7 r 

sm ft = \ — and 

1 ' V 2 124 82 



1 . . 1 l + y/sin^f-sin^) 4 . ^ 



sin(6»i) = \ - + -i — 

1 ' \2 37446 



785323439 

and 



I y + v/3 J -I + 28sin(6»i) 2 - 28sin(6»i) 4 

1 > sm(6 2 ) > -— = \ ; 77~\2~T • Tk .4 and 1 > sin (#3 

3y2 \ 1 — sin(6q) + sm(#i) 



> 



V3 



— = < sin(^i) < \ - + 

V ' \ 2 37446 



785323439 

and 



79 
100 

1 



<sin(0 2 )<-^ 



9 + V3 



1 - sin(6' 1 ) 2 + sm(6 



i 2 - 28sin(0i) 4 

4 



-l + 28sin(6>i) 

and 



1 > sm(6 3 ) > -— \ 



3 + V3 y/-l + 4sin(6» 2 ) 2 - 4sin(6l 2 ) 4 + 4sin(0i) 2 sin^) 4 - 4sin(6»i) 4 sin(6> 2 ) 4 



1 — sin^) + sin(# 2 ) — sinl 



9i) 2 sin(0 2 ) 4 + sin(0 1 ) 4 sin(0 2 ) 



(58) 



(59) 



(60) 
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1 > sin 




V32789757 
and 



79 

1 > sin(0 2 ) > Y?j7j and 1 > sin lf:; 



(61) 



V32789757 
12482 




> sm(6>i 



785323439 



37446 



and 1 > sin( 



79 

h) > — and 1 > sin( 



^3 
2 



> ( 62 ) 



1 > sin(0 2 ) > - 



V2 



< sin(6»i) < 



785323439 

and 



37446 



i 

3^2 \ 



VW-1 + 28sin(6»i) 2 - 28sin(6»i) 4 ^ 

* -, and 1 > sin(# 3 ) > -— . 

1 - sin(0i) 2 + sin(^) 4 2 



(63) 



Therefore we could generate a pd which could be entangled but would no longer be a pure state. Unfortunately, the 
fact that pd is no longer a pure state also means that we would have to look at the most general U € SU (4) acting on 
Pd in order to determine which successive unitary operations UpdU^ would produce entanglement. This is a rather 
lengthy and complicated calculation and is beyond the scope of this section. 

What we can do though is make an educated guess as to the volume of entangled two qubit mixed states by noticing 
that it is a product of the volume of the 3-dimensional symplex of eigenvalues of pd (with appropriate ranges) and 
the volume of the flag manifold S[/(4)/(/(l) S(/ (^i x U(l)su(3) x U(l) su ^y 10 These volumes can be calculated by 
using the mixed state product measure given in [3j with the necessary ranges for the eigenvalues (given above) and 
the covering ranges for ST/ (4) (given in Using this material we arrive at 



dV Ems =a s A s l - 1 A s 2 - 1 A s 3 - 1 (l-^2Ai) s - 1 dAi...dA 3 x d 



517(4) 



K U(1) SU {2) X U(l) su{3) X U(1) SU (4) 

A general (and rather naive) evaluation of this measure for our situation yields 



V En 



(a L > - a v s ) (b L s - b v s ) (c L s - c v s ) (d L s - d v s ) vr 6 tt 
a s a x — = uj 



12 



12' 



da\2 ■ ■ ■ da\. (64) 



(65) 



where {ajj, ol},{6(7, 6l},{c[/, cl} and {o?[/,<7l} are the squared values of the above maximal and minimal ranges 
({djj, c?l} comes from A 4 = 1 — 53i=i ^»)- Since the symplex measure is assumed to be non-zero, and using the work 
contained in Q we can hypothesize that u> has the following bounds (dependent on the value of s and recalling that 
a s > and s > 0) 



< (a L s ~ a v s ) (b L s - b v s ) (c L s - c v s ) (d L s - d v s ) < 4- 4s (-l + 4 s ) 



(66) 



Notice that the right side of the above inequality approaches when both s — > oo and s — > 0, thus we can conclude 
that the numerator of (|65l) will be < 1 and therefore, the value of w (be it either > 1 or < 1) will be completely 
dependent on the explicit choice of the value of s. 

The important point to recognize here is not the symplex calculations but rather the flag manifold volume. Through 
the Euler parameterization of SU(N) and U(N) given in 0, Q we have been able to generate the appropriate 
representation of the "truncated" Haar measure which is crucial to any mixed state volume calculation. It is this 
factor which is not "user dependent"; i. e. dependent on the initial distribution chosen for the (N — l)-dimensional 
symplex, and therefore not completely subject to disagreements between researchers studying entanglement. 11 



10 Explained in detail in 0. 

11 Although disagreements in numerical values are found, they are mostly due to variations in the ranges of the N(N — 1) parameters 
which define the measure (see for example |14| and references within). 
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III. QUBIT/QUTRIT ENTANGLEMENT 
A. Pure State Entanglement 



By following the same procedure as was done in the two qubit case, we can derive the manifold of operations that 
produce entanglement of an initial qubit / qutrit pure state 



Pd 



/l 0\ 









\o 0/ 



(67) 



via a U G SU(6) given in 



U = U (ai, a 2 , a 3 , <xi, a 5 ,a 6 ,a 7 , a 8 , ag, aio, an, 0, ai 3 , 0, a%5, 0, an, 
x 0, ai 9 , 0, a 2 i, 0, a 23 , 0, a 25 , 0, a 27 , 0, a 2 g, 0, a 3 i, a 32 , a 33 , a 3 4, ^35) 

£j _ giA 3 Qi e iA2a2giA3a3 e iA5a4giA3Q5 e iAioa6 e iA3Q7 e iAi7a 8e iA 3 a9 e iA26aio 
„ „jA3Cm iA3C«i3 iAsais iA3Cti7 iA3ai9 iA3U2l „iA3Q;23 

AC C C t C C C 

x giA 3 Q;25 giA 3 Q 2 7 giA 3 a29 giA 3 a 3 i gjAi 5 a 3 3 giA 2 4Q34 gjA 35 a35 



(68) 



by taking the partial transpose of U pjU^ and evaluating the corresponding eigenvalues through the Peres-Horodccki 
criterion. 12 Doing this work, which did not yield well to simplification, generates the following hypothesized manifold 



SU{6) 



CP 



Me ^< 6 > U(5) x U(l) SU (2U U(1) SU(2) J 



(69) 



where k 



-2(2ai + a 3 + (17 + ag) is the cumulative phase. 13 The volume of the manifold is 

SU(6) \ V r 



Vj 



V 



U(5) x U(1) su(2)k 



10* Vr 



(7(1) 



5! * 20tt 4 

7T 



2400 



(70) 



If we make the same argument for the volume of this manifold to be equivalent to the volume of entangled qubit / qutrit 
pure states, as we did in the two qubit case, then, in this case, our values are within the ranges specified in but 
since we do not include the possibility of bound entangled states (states which are entangled, but which have positive 
partial transposes) we must concede that our value given in equation l|70|) is probably too small. 

B. Mixed State Entanglement 

Again, as in the two qubit case, one can see that for a general pd for a qubit/qutrit system 

n(6 , 3 ) 2 sin(6l4) 2 sin(6'5) 2 , 



p d = diag{sin((9i) 2 sin(6> 2 ) 2 sin(i, 3 , 0111^4; c 
s(0 2 ) 2 sin(0 3 ) 2 sin(04) 2 sin(0 5 ) 2 , 



cos( 



:) 2 sin(6> 4 ) 2 sin(# 5 ) 2 , cos(#i) 2 sin(0 2 ) 2 sin(# 3 ) 2 sin^) 2 sin(6 , 5 ) 2 , 
t) 2 sin(0 5 ) 2 , cos(# 3 ) 2 sin(^4) 2 sin(0 5 ) 2 , cos(^4) 2 sin(6> 5 ) 2 , cos^) 2 } 



(71) 



We are forced to generate the six eigenvalue equations rather than evaluate the constant term from the characteristic polynomial because 
in this case we have two eigenvalues equal to zero thus negating the constant term's effectiveness. 

After taking into account the degeneracy in the eigenvalues, full simplification of the resulting 4th order characteristic polynomial was 
not possible without making certain numerical assumptions. Thus, the actual representation of k as a function of Qi, Q3, c?7 and og was 
not possible; thus the equivalence. 
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when 6i 7^ tt/2, (i = 1, . . . , 5) the generalization of equation (|13|) for the qubit/qutrit case 



5 

where 



T ? V] > \ (72) 



Tr[p 2 ] ^TrKC/p.f/t) 2 ] = Tr[p 2 d ] 

= 1 - 2sin((9 5 ) 2 (1 + sin(0 5 ) 2 - sin(6» 4 ) 2 sin(6> 5 ) 2 + sm(6» 4 ) 4 sin(6» 5 ) 2 

— sin(# 3 ) 2 sin(^4) 4 sin(6> 5 ) 2 + sm{Q s ) A sin(6*4) 4 sin^) 2 — sin(0 2 ) 2 sin(# 3 ) 4 sin(#4) 4 sin^s) 2 

+ sin(# 2 ) 4 sin(# 3 ) 4 sin(6*4) 4 sin^) 2 — sin(6*i) 2 sin^) 4 sin(# 3 ) 4 sin(#4) 4 sin(# 5 ) 2 

+ sin(6» 1 ) 4 sin(6» 2 ) 4 sin((9 3 ) 4 sin(6> 4 ) 4 sin((9 5 ) 2 ) (73) 

could still be satisfied. In this case then we would have for the entangled mixed state product measure (under 
appropriate ranges for Aj and cti) 

5 

dV Ems =a s A s 1 ~ 1 ■ ■ ■ A. s 5 -\l - A l ) s ~ 1 dAi ...dA 5 

i=l 

J SU(6) \ , 

x a — —- ^- t— r da^o . . . dot\ . (74) 

\U(1)SU(2) X U{1) SU(3) X U{1) SV (4) X U{l)g U{5) X [7(1) 

SU(6) / 

Another general (and again, rather naive) evaluation of this measure for our situation yields 

(a L s - a v s ) (b L s - b v s ) {c L s - c v s ) (d L s - d v s ) (e L s - e v s ) (f L s - f v s ) tt 15 



V Ems = a s 



34560 



TT 15 



34560' 



(75) 



where {ajj, a^} through {/[/, fjj} are the squared values of the maximal and minimal ranges of A, that satisfy l|72f) . 
As before, the symplex measure must be assumed to be non-zero, therefore using the work contained in Q we can 
hypothesize that uj has the following bounds (dependent on the value of s and recalling that a s > and s > 0) 

< (a L s - a v s ) (b L s - b v s ) (c L s - c v s ) (d L s - d v s ) (e L s - e v s ) (f L s - /V) < 6~ es (-l + 6 s ). (76) 

Again we notice that the right side of the above inequality approaches when both s — > oo and s — > 0, we can 
again conclude that the numerator of 175|) will be < 1 and therefore, the value of u> (be it either > 1 or < 1) will be 
completely dependent on the explicit choice of the value of s. Also, as before in the two qubit case, the important 
point to recognize here is not the symplex calculations but rather the flag manifold volume. 



IV. CONCLUSIONS 



In this paper we have applied our SU{N) and U{N) parameterizations to the two qubit and qubit/qutrit system 
in order to explicit calculate the manifold of operations which entangle two qubit and qubit/qutrit pure states. We 
have also been able to give the volume of this manifold, as well as the hypothesized volume for the set of all entangled 
two qubit and qubit/qutrit pure and mixed states. In the pure state case, the values were within the ranges given 
by 1^,111 but in the qubit/qutrit case, because we did not take into account the possibility of bound entangled states 
(which do not appear in the two qubit case) our volume is most likely smaller than the actual volume for the set of 
all entangled qubit/qutrit pure states. 

Work is continuing on the mixed state situation; explicitly in calculating the volume of the mixed state manifold 
without having to know the exact probability distribution on the [N — l)-dimensional symplex. Extensions of the 
pure state work to two qutrit systems is also ongoing. 
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APPENDIX A: SU(i) LIE ALGEBRA 



5] 



From Q] we know that the Gell-Mann type basis for the Lie algebra of SU(4:) is given by the following set of matrices 



Ai = 



A 4 = 



A 7 = 



Aio — 



U3 




As = 




A 3 = 



A fi = 



Al2 — 



A IF 




(Al) 



Using these matrices one can then generate the various group operations given in section [n] Similarly, in f| one can 
see how to construct the N 2 — 1 elements of the SU(N) Lie algebra necessary for general SU(N) group operations. 
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